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Abstract 


This  paper  studies  the  dynamics  of  coupled  planar  rigid  bodies,  concentrating  on  the 
case  of  two  or  three  bodies  coupled  with  a  hinge  joint.  The  Hamiltonian  structure  is  non- 
canonical  and  is  obtained  using  the  methods  of  reduction,  starting  from  canonical  brackets 
on  the  cotangent  bundle  of  the  configuration  space  in  material  representation.  The  dynamics 
on  the  reduced  space  for  two  bodies  occurs  on  cylinders  in  stability  of  the  equilibria  is 
studied  using  the  Energy  -  Casimir  method  and  is  confirmed  numerically.  The  phase  space  of 
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the  two  bodies  contains  a  homoclinic  orbit  which  produces  chaotic  solutions  when  the  system 
is  perturbed  by  a  third  body.  This  and  a  study  of  periodic  orbits  are  discussed  jn  part  D.  The 
number  and  stability  of  equilibria  and  their  bifurcations  for  three  bodies  as  system  parameters 
are  varied  are  studied  here;  in  particular,  it  is  found  that  there  are  always  4  or  6  equilibria. 
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§1.  Introduction 

The  techniques  of  reduction  of  Hamiltonian  systems  with  symmetry  and  the  attendant 
Energy  -  Casimir  method  have  proved  to  be  useful  in  a  wide  variety  of  problems,  including  fluid 
and  plasma  stability  (Holm,  Marsden,  Ratiu  and  Weinstein  [l985]),  rigid  body  dynamics  with 
attachments  and  internal  rotors  (Holmes  and  Marsden  [1983],  Koiller  [1985],  Krishnaprasad 
[1985],  Krishnaprasad  and  Marsden  [1987]),  and  bifurcations  of  liquid  drops  (Lewis,  Marsden 
and  Hatiu  [1986a, b]).  In  this  paper  we  shall  apply  these  techniques  to  the  case  of  planar  rigid 
bodies  coupled  by  a  hinge  joint.  Many  of  the  results  for  the  two  and  three  bodies  generalize 
to  multibody  structures  and  other  modifications,  such  as  the  inclusion  of  hinge  torques.  In 
subsequent  papers  we  shall  be  studying  this  as  well  as  the  problem  of  coupled  three  dimensional 
rigid  bodies  (e.g.,  with  a  ball-in-socket  or  hinge  joint).  We  also  expect  that  the  non-canonical 
Hamiltonian  methods  that  are  useful  here  will  also  be  useful  in  related  problems  of  control 
(see  Van  der  Schaft  [1984]  and  Sanchez  de  Alvarez  [1986]). 

The  reduction  technique  used  here  goes  back  to  Arnold  [1966],  Meyer  [1973],  and  Marsden 
and  Weinstein  [1974],  amongst  others.  It  involves  starting  with  a  Poisson  manifold  P  and  a  Lie 
group  C  acting  on  P  by  canonical  transformations.  The  reduced  phase  space  P/G  (assume  it 
has  no  singularities)  has  a  natural  Poisson  structure  whose  symplectic  leaves  are  the  Marsden- 
Weinstein-Meyer  spaces  J~^(fi]/Gtt  «  J~^{0)/G  where  (t  £  g*,  the  dual  of  the  Lie  algebra 
of  G,  J  :  P  — »  y*  is  an  equivariant  momentum  map  for  the  action  of  G  on  P,  G^  is  the 
isotropy  group  of  p  (relative  to  the  coadjoint  action)  and  0  is  the  coadjoint  orbit  through  p. 
If  P  =  T*G  and  G  acts  by  left  translations,  then  P/G  is  identifiable  with  y*  equipped  with 
the  (•)  Lie-Poisson  bracket: 

{f.FKo)  =  -<»*.[^.^|>  (11) 

The  symplectic  leaves  in  this  case  are  just  the  coadjoint  orbits.  For  G  =  SO (3)  we  get  the 
(Pauli-Martin)  bracket  for  rigid  body  dynamics: 

{P,P}(/)  «  -t.(VFxVB)  (1.2) 

Here  t  €  bo(3)*  is  identified  with  a  vector  in  and  represents  the  angular  momentum  of  the 
rigid  body  in  a  body-fixed  frame.  If  1  is  the  moment  of  inertia  tensor  so  ^  =  I  w  where  w  is 


S 


tbe  body  angular  velocity,  then  Euler^a  equations: 


^  txiif  .  (1.3) 

are  equivalent  to  Hamilton’s  equations 

^  “  {F,B)  (1.4) 

where  H(/)  =  |<£,ii;>= 

Notice  that  (1.2)  is  a  non-canonical  bracket;  i.e.  the  usual  (9,p)  Poisson  bracket  formalism 
has  disappeared  through  the  reduction  process.  One  of  our  first  goals  in  the  paper  will  be  to 
develop  a  similar  bracket  for  the  dynamics  of  two  coupled  planar  rigid  bodies.  We  start  with 
the  canonical  bracket  on  the  cotangent  bundle  of  configuration  space  just  as  one  starts  with 
T*SO(S)  (parametrized  by  Euler  angles  (d,^,^)  and  their  conjugate  momenta  (p#iP^,PV')) 
in  rigid  body  dynamics. 

When  these  procedures  are  carried  out  for  coupled  rigid  body  dynamics  (§§2*4)  we  find 
that  concepts  akin  to  the  ’augmented  body’  (cf.  Wittenburg  [1977])  come  out  in  a  natural 
way.  The  reduced  Poisson  structure  obtained  is  a  Poisson  structure  in  R’  (not  of  Lie>Poisson 
type,  however)  whose  symplectic  leaves  are  cylinders.  The  reduced  dynamics  on  one  of  these 
cylinders  for  specific  rigid  body  parameters*  is  shown  in  Figure  1.  Here,  fii,  m  are  closely 
related  to  the  angular  momenta  of  the  two  bodies  and  $  is  the  joint  angle. 


*  Tbe  parameters  chosen,  in  the  notation  of  i|2'4  are  Ji  »  105.S5,  />  s  70,  c  ^  5S.55,  and 
Ml  +  Ma  *  *0- 
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Being  two  dimensional  and  Hamiltonian,  the  flow  on  the  cylinder  is  completely  integrable. 
Notice  that  there  are  two  equilibria,  one  a  saddle  and  one  a  stable  point.  This  is.conflrmed  by 
a  linearized  analysis  for  the  saddle  point  and  an  Energy  •  Casimir  analysis  for  the  stable  point 
(see  Holm  et  al  [1985]).  The  stable  point  corresponds  to  the  two  bodies  uniformly  rotating 
in  an  extended  position,  while  the  saddle  point  corresponds  to  uniform  rotation  in  a  folded 
position  (figure  2) 


body  1 


body  2 


(a)  stable  equilibrium 


body  1 


(b)  unstable  equilibrium 


Figure  2 


There  are,  of  course,  corresponding  equilibria  for  oppositely  oriented  rotational  motions. 

Notice  from  Figure  1  that  there  are  two  homoclinic  orbits  from  the  unstable  equilibrium 
back  to  itself.  Thus,  one  can  expect  that  when  for  example,  an  additional  third  body  is 
attached  nearly  at  the  center  of  mass  of  body  2  or  the  system  is  forced  (for  instance  by  joint 
torques),  there  will  be  a  splitting  of  these  bomoclinic  orbits  resulting  in  chaotic  dynamics. 
One  way  to  proceed  with  an  analysis  of  this  sort  is  via  the  Melnikov  method  (see  Holmes 
and  Marsden  [1982,1983]  and  Guckenheimer  and  Holmes  [1983]).  This  analysis,  together  with 
more  information  on  instability  and  periodic  orbits  will  be  given  in  part  H  of  this  paper. 

Another  benefit  of  doing  the  analysis  systematically  using  the  reduction  procedure  is  that 
the  generalization  to  multibody  problems  and  three  dimensional  motion  can  be  done  using 
similar  ideas.  We  discuss  the  planar  multibody  case  in  $6  and  the  three  dimensional  case  in 
another  publication. 


We  now  summarize  one  of  the  results  of  the  present  work;  namely  we  display  the  Hamilto¬ 
nian  form  for  the  dynamic  equations.  The  details  of  the  derivation  of  this  structure  are  given 
in  §§2-4.  Refer  to  Figure  3  and  define  the  following  quantities: 


Figure  3 


di  =  distance  from  the  binge  to  the  center  of  mass  of  body  «  =  1,2 
u)i  =  angular  velocity  of  body  t  =  1,2 
9  =  joint  angle  from  body  1  to  body  2 
A(0)  =  did]  cos 9 
nti  =  mass  of  body  t  =  1,2 
e  =  mim3/(mi  +  ma)  =  reduced  mass 
Ji  =  moment  of  inertia  of  body  i  about  its  center  of  mass 

/i  =  A  *4-  «dj  ;  /j  =  /j  +  cdj  s=  augumented  moments  of  inertia. 


and  7 


c  V 

(/,/,  -  c»A»j 


The  dynamics  of  the  system  is  described  by  the  following  Euler-Lagrange  equations  for 
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#,  «!,  wa: 


9  =  —  ux 

wi  =  -7(/aw|  + 

^2  =  +  <Aa;|) 

For  tbe  Hamiltonian  structure  it  is  convenient  to  introduce  the  momenta 


/ii  =  ii^x  +  cXu>2 

i.e. 

H2  =  /jwj  +  cAwi 


(1.5) 


(1.6) 


(this  b  done  via  the  Legendre  transform  in  §4.)  The  evolution  equations  for  Hi  are  obtained 
by  solving  (1.6)  for  Ux,  ux  and  substituting  into  (1.5).  The  Hamiltonian  is 

(;;;*)  (w.) 

i.e. 

which  b  the  total  Icinetic  energy  for  the  two  bodies.  Tbe  Poisson  structure  on  the  (^,/xi«M2)* 
space  (called  P  in  §3)  b 

{F.H}  =  {F.H),  -  {F.fl),  (1.8) 

where 

irm  - 

I/-,///.  _  . 

The  evolution  equations  (1.5)  are  then  equivalent  to  Hamilton’s  equations  F  =  {F,H}. 
Casimirs  for  tbe  bracket  (1.8)  are  readily  checked  to  be 

C  =  ♦(m,  +  Ms)  (1-9) 

for  ft  any  function  of  one  variable;  i.e.  {F^C)  s  0  for  any  F.  One  can  also  verify  directly  from 
(1.5)  that,  correspondingly,  =  0  where  m  =  Mi  +  Ms  t**®  system  angular  momentum. 
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The  symplectic  leaves  of  (1.8)  are  described  by  the  variables  v  —  {H2  /^i)/2,8  which 

parametrize  the  cylinder  shown  in  Figure  1.  The  bracket  in  terms  of  {9,v)  is  the  canonical 


one  on  T*S^: 


{F,H}  = 


dF  an 

89  du 


80  8u 


(1.10) 


As  we  shall  see,  this  canonical  structure  on  T'S*  is  consistent  with  the  Satzer-Marsden- 
Kummer  cotangent  bundle  reduction  theorem  (Abraham  and  Marsden  [1978],  Kummer  [1981]). 


§2.  Kinematlcal  Preliminaries  (for  two  coupled  planar  rigid  bodies) 

In  this  section  we  set  up  the  phase  space  for  the  dynamics  of  our  problem.  Refer  to  Figure 
4  and  dehne  the  following  quantities: 


"Ty 


center  of  mass 
of  the  system 


reference 

point 


Figure  4 


di3  =  the  vector  from  the  center  of  mass  of  body  1  to  the  binge  point  in  a  reference  config¬ 
uration  (fixed). 

dsi  =  the  vector  from  the  center  of  mass  of  body  2  to  the  binge  point  in  a  reference  config¬ 
uration  (fixed). 

i?(tf^)  s=  cos°^^*)  ~  rotation  through  angle  9i  giving  the  current  orientation  of 

body  I  (written  as  a  matrix  relative  to  the  fixed  standard  inertial  frame), 
r*  =  current  position  of  the  center  of  mass  of  body  t 
r  =  current  position  of  the  system  center  of  mass 
T®  =  the  vector  from  the  system  center  of  mass  to  the  center  of  maiss  of  body  i. 

#  =  02-^1  =  joint  angle 
JJ(#)  =  joint  rotation  =  J2(tfa)  •  i?(— #j) 

The  basic  configuration  space  we  start  with  is  Q,  the  subset  of  5.E(2)  x  5E(2)  (two 
copies  of  the  special  Euclidean  group  of  the  plane)  consisting  of  pairs  ((i2(tfj),ri),  (E(02)<r3)) 


0 


satisfying  the  hinge  constraint 


Fj  =  Ti  +  Ji(fii)di2  —  ^(^2)^21  (2-1) 

Notice  that  Q  is  of  dimension  4  and  is  parametrized  by  ^1,^2  and,  say  Fj;  i.e.  <3  «  5*  xS*  xR*. 
We  form  the  velocity  phase  space  TQ  and  momentum  phase  space  T*Q. 

The  Lagrangian  on  TQ  is  just  the  kinetic  energy  (relative  to  the  inertial  frame)  given  by 
summing  the  kinetic  energies  of  each  body.  For  convenience,  we  recall  how  this  proceeds:  let 
Xi  denote  a  position  vector  in  body  1  relative  to  the  center  of  mass  of  body  1,  and  let  Pi(Xj) 
denote  the  mass  density  of  body  1.  Then  the  current  position  of  the  point  with  material  label 


Xi  is 


Thus 


Xi  =  .R{tfi)Xj  +  Fj 


xi  = 


and  so  the  kinetic  energy  of  body  1  is 


=  II  X.  f  i'X. 

~lj  +  *1  >  <f’Xi 

=  l/(n(X,)|<  fiXi.AX,  >  +2<AX„r,>  +  llTill’li’: 

<flX.,AX,  >=  (r(flX.(KX.)’')  =  <r(j?xrx,«’') 


f„{X,)  <  JbCi.r,  >  i’X,  =  <  ,.CX.)X.<l”X.,r.  >  =  0  (2.5) 

since  Xi  is  the  vector  relative  to  the  center  of  mass  of  body  1.  Substituting  (2.4)  and  (2.5) 
into  (2.3)  and  defining  the  matrix 

V  ~  j  l>(X,)X,XT<l'X,  (2 «) 


Ki  = 


(2.7) 


|(r(K(«Ol'«(».)’')  +  Im,  II  r.  II’ 

with  a  similar  expression  for  we  let 

L:  TQ^VL  be  L  =  Ki  +  (2.8) 

The  equations  of  motion  then  are  the  Euler-Lagrange  equations  for  this  L  on  TQ.  Equivalently, 
they  are  Hamilton's  equations  for  the  corresponding  Hamiltonian. 

For  later  convenience,  we  shall  rewrite  the  energy  (2.8)  in  terms  of  =  6i,  wj  =  $2iTi 
and  rf.  To  do  this,  note  that  by  definition, 

mr  =  miTi  +  m2T2  (2.9) 

where  m  =  mj  +  mj,  and  so,  as  ri  =  r  +  r®, 

0  =  miT®  +  mjrf  (210) 

and,  subtracting  r  from  both  sides  of  (2.1), 

r®  =  r®  +  -R(tfi)d:2  -  RiWii  (2.11) 

From  (2.10)  and  (2.11)  we  find 

rS  =  —  (i2(tfi)dxj  -  i?(«2)d2a)  (2.12a) 

tn 

and 

r®  =  !!^(i2(tf,)dj2  -  i?(^2)d2i)  (2.126) 

tn 

Now  we  substitute 


into  (2.8)  to  give 

2 

+  ilm.(l|r  +  iff)  +  m,  Ilf  +  rill*) 


+  jlm.dl  t  +  f{  1|»)  +  m,  11  f  +  fi  ll’l  (2.14) 

2 

<  r,r5  >=  0  since  mii^  +  mjfS  =  0  from  (2.10).  Thus  (2.14) 


But  mi  <  r.fj  >  +  ma  <  r,r5  >=  0  since  miiS  +  mafg  =  0  f 
simplifies  to 

i=  i<r(K(f.)i'fi(«.)’'  +  6{e2)l’R(hf) 
2 


^  +  |m.||fi||*  +im.||ff||* 


where  p  =  m 


2m  2  * 

T  II  is  the  magnitude  of  the  system  momentum. 


Now  write 

£(cosei  -Bintfi\ 
dt\sln$i  cos«i  ; 


1  that  (2.12)  gives 


dt  \  sin  $i  cos  9i  ) 

=  ('"V;  ==  ’o’)  = 

es 


fO  =  ^(i2(tfl)ildj2  -  J2(^2)w2d2l) 


m 

fO=  ^(iZ(tfi)widia  -  R{h)^2^2i) 


Substituting  (2.17)  and  (2.16)  into  (2.15)  gives 

1  .  «v*.  ^ 


Finally  we  note  that 


L=  5tr(wil*«ir)  +  W2l*wf)) 

2 

+  i  II  6,d„  -  * (»2  -  II’ 

2m  m 


(2.15: 


(2.16 


(2.17a 


(2.176 


wbere 


/i  =  j  p{Xi,Yx){Xl  +  Y,^)dXxdYi 

is  the  usual  moment  of  inertia  of  body  one  about  its  center  of  mass.  One  similarly  derives 
(2.192]  is  which  “1"  is  replaced  by  “2".  The  last  term  in  (2.18)  is  manipulated  as  follows: 

II ‘*'1^12  “  i2(^)<^2d2l  II*  =  11  Wldj2  11*  —  2  <  Widi2,i?(^)u>2d2l  > 

+  11  wjdai  II* 

=  u^di  +  u\d2  —  2  <  iuxdi2t(^2'R(^)d2i  > 

=  C^*dj  +  b>2^2  ~  2tUiU»2  <  di2» /Z(^)d21  >  (2.20) 

Substituting  (2.19}),  (2.192)  and  (2.20)  into  (2.18)  gives 

Xr  =  -  ^(w*/i  +  ‘*'1^2)  +  2cJiU2eX{9)  +  (^-^l) 

where 

A(^)  =  “  <  di2, •R(^)d2i  >  =  — |di2  *  ^21  cos  --  (dx2  ^  d2i)  *  ^sin9]  (2.22) 

Remarks:  1.  If  dx2  d2i  are  parallel  (i.e.  the  reference  configuration  is  chosen  with  dx2 

and  d2x  aligned),  then  (2.22)  gives  A(0)  s  (iid2COs9,  as  in  $1. 

2.  /x,  I2  are  the  moments  of  inertia  of  "augmented'’  bodies  as  defined  in  §1;  for  example 
/x  is  the  moment  of  inertia  of  body  1  augmented  by  putting  a  mass  c  at  the  hinge  point. 


§3.  Reduction  to  the  Center  of  Mass  FVame 

In  this  section  we  reduce  the  dynamics  by  the  action  of  the  translation  group  R’.  This 
group  acts  on  the  original  configuration  space  Q  by 

v.((i?(tfx),Pi),(R{tf2).'2))  -  ((/2(«i),ri+v).(R{tf2),ra+v))  (3.1) 

This  is  well-defined  since  the  hinge  constraint  (2.1)  is  preserved  by  this  action.  The  induced 
momentum  map  on  TQ  is  calculated  by  the  standard  formula 

=  IlfaW  (’-S") 


or  on  T*Q  by 

Jt  =  Pi£e(«) 

where  {^is  the  infinitesimal  generator  of  the  action  on  Q.  (See  Abraham  and  Marsden  [1978]). 
To  implement  (3.2)  we  parametrize  Q  by  9t,  $z  and  r  with  Pj  and  determined  by  (2.12) 
and  (2.13).  From  (2.15)  we  see  that  the  momentum  conjugate  to  r  is 

p  =  ^  =  mf  (3.3) 


and  so  (3.2)  gives 


=  <  P*(  >»  C  G 


(3.4) 


Thus  <7  =  p  is  conserved  since  H  is  cyclic  in  r  and  so  H  is  translation  invariant.  The  corre¬ 
sponding  reduced  space  is  obtained  by  fixing  P  =  Po  letting 


Pro  -  J-‘(Po)/R’ 


(See  Abraham  and  Marsden  (1978, Ch.4]).  But  is  clearly  isomorphic  to  r*(5*  x  5*)  i.e. 
to  the  space  of  #j,  #3  and  their  conjugate  momenta.  The  reduced  Hamiltonian  is  simply  the 
Hamiltonian  corresponding  to  (2.21)  with  p  regarded  as  a  constant. 

Note  that  in  this  case  the  reduced  symplectic  manifold  is  a  cotangent  bundle,  in  agreement 
with  the  cotangent  bundle  reduction  theorem  (Abraham  and  Marsden  [1978],  Rummer  [1981]). 


can  also  check  this  directly 


The  reduced  phase  space  has  the  canonical  symplectic  form;  one 

here. 

In  (2.21)  we  can  adjust  L  by  a  constant  and  thus  assume  p  =  0;  this  obviously  does  not 
affect  the  equations  of  motion. 


hti  us  observe  that  the  reduced  system  is  given  by  geodesic  flow  on  S'  x  S'  since  (2.21) 
is  quadratic  in  the  velocities.  Indeed  the  metric  tensor  is  just  the  matrix  J  given  by  (1.6),  so 
the  conjugate  momenta  are  pii  P2  given  by  (1*6). 


We  remark,  finally,  that  the  reduction  to  center  of  mass  coordinates  here  is  somewhat 
simpler  and  more  symmetric  than  the  Jacobi-Haretu  reduction  to  center  of  mass  coordinates 
for  n  point  masses.  (Just  taking  the  positions  relative  to  the  center  of  mass  does  not  achieve 
this  since  this  does  not  reduce  the  dimension  at  all!)  What  is  different  here  is  that  the  two 
bodies  are  hinged,  and  so  by  (2.12),  and  ra®  »re  determined  by  the  other  data. 
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§4.  Reduction  by  Rotationt 

To  complete  the  reduction,  we  reduce  by  the  diagonal  action  of  S*  on  th^  configuration 
space  X  that  was  obtained  in  §3.  The  momentum  map  for  this  action  is  obviously  given 

by 

=  Ml  +  /*a  (4'1) 

For  purposes  of  later  stability  calculations,  we  shall  find  it  convenient  to  form  the  Poisson 
reduced  space 

P  ;=  r*(S*  X  S') /S'  (4.2) 

whose  symplectic  leaves  are  the  reduced  symplectic  manifolds 

=  J^'M/s'cP 

We  coordinatize  P  hy  $  =  9^  -  Ml  Mai  topologically,  P  =  S*  x  R*.  The  Pois¬ 
son  structure  on  P  is  computed  in  the  standard  way:  take  two  functions  P(^,Mi»Ma)  &Dd 
P(^,Mi,Ma).  Regard  them  as  functions  of  #i,#3,MiiMa  by  substituting  9  -  9^  -  9i  and 
compute  the  canonical  bracket.  It  is  clear  that  the  asserted  bracket  (1.8)  Is  what  results.  The 
Casimirs  on  P  are  obtained  by  composing  J  with  Casimirs  on  the  dual  of  the  Lie  algebra  of 
S*;  i.e.  with  arbitrary  functions  of  one  variable;  thus  (1.9)  results.  This  can  of  course  be 
checked  directly. 

If  we  parametrize  P^  by  9  and  v  =  then  the  Poisson  bracket  on  P^  becomes  the 

canonical  one.  This,  again,  b  consistent  with  the  cotangent  bundle  reduction  theorem  which 
asserts  in  this  case  that  the  reduction  of  r*(S*  x  S*)  by  S*  b  symplectically  diffeomorphic 
to  T*{{S'  X  5*)/5*)  S?  T*S'.  There  are  no  'magnetic’  terms  since  the  reduced  configuration 
space  s'  b  one  dimensional,  and  hence  has  no  non-sero  two  forms. 

The  realization  of  P^  as  7*5*  b  not  unique.  For  example  we  can  parametrize  P^  by 
(^a.Ma)  or  by  (^i,Mi)t  «*ch  of  which  also  gives  the  canonical  bracket.  (In  the  general  theory 
there  can  be  more  than  one  one-form  “an*  by  which  one  embeds  P^  into  r*5*,  as  well  as 
more  than  one  way  to  identify  (S*  x  S')/S'  K  S'.  The  three  listed  above  correspond  to 
three  such  choices  of  tti,)* 

Remark.  The  reduced  bracket  on  T*{S'  xS')/S'  can  also  be  obtained  from  the  general  formula 
for  the  bracket  on  (P  x  T*C)/G  S£  P  x  g*  found  in  Krishnaprasad  and  Marsden  (1988];  it 
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produces  one  of  the  variants  above,  depending  on  whether  we  take  G  to  be  parametrized  by 
#x  or  #2,  or  ^3  —  ^i. 

The  reduced  Hamiltonian  on  P  is  just  (l.7b)  regarded  as  a  function  of  MiiMs  We 

therefore  know  that  the  Euler- Lagrange  equations  (1.5)  are  equivalent  to  P  =  {P,P}  for  the 
reduced  bracket  (1.8). 

We  can  also  obtain  a  Hamiltonian  system  on  the  leaves,  parametrized  by  say  {9,u).  We 
simply  take  (1.7b),  namely 

T](«) 

where  A  =  /1/2  -  <’A*,  and  substitute  (ii  =  f  -  t'*  Ma  =  u  +  f  producing 

+  +  /»-*<>))  (40 

The  presence  of  the  linear  term  in  1/  can  be  eliminated  by  completion  of  squares:  it  b  not 
there  in  the  general  theory  (Abraham  and  Marsden  [1978],  Smale  [1970])  because  reduced 
coordinates  adapted  to  the  metric  of  the  kinetic  energy  are  used;  these  are  produced  by  the 
completion  of  squares.  Notice  that  the  Hamiltonian  now  is  the  form  of  kinetic  plus  potential 
energy  but  that  the  metric  now  on  is  8-dependent  and,  unless  di  or  dz  vanishes,  it  is  a 
non-trivial  dependence.  The  potential  piece  is  usually  referred  to  as  the  amended  potential 

We  summarize  : 

Theorem  1  :  The  reduced  phase  space  for  two  coupled  planar  rigid  bodies  is  the  three 
dimensional  Poisson  manifold  P  =  S*  x  R  with  the  bracket  (l-B);  sympleetic  leaves  are  the 
cylinders  with  canonical  variables  (8,1').  Casimirs  are  given  by  (1-8). 

The  reduced  dynamics  are  given  by  P  {P,  B)  or  equivalently, 


vihtre  E  is  given  by  (1.7b).  The  equivalent  dynamics  on  the  leaves  is  given  by 


§5.  Equilibria  and  Stability  by  tbe  Energy  >  Casiinir  Method 


We  now  use  Arnold’s  energy-CasImir  method,  as  is  summarized  in  Holm  et.al.  [1985]  and 
Krishnaprasad  and  Marsden  [1986]  to  determine  the  equilibrium  points  and  their  stability.  An 
equivalent  alternative  to  this  method  is  to  look  for  critical  points  of  H  given  by  (4.4)  in  ($,  i/) 
space  and  test  for  definiteness  at  these  equilibria. 

To  search  for  equilibria  we  look  directly  at  Hamilton’s  equations  on  P.  Using  the  bracket 
(1.8)  and  F  =  {F,E},  we  obtain  equations  (4.5),  where  B  is  given  by  (l.7b).  The  conditions 
fii  pt2  =  0  become 


dB 

as 


0; 


i.e. 


=  0 


Clearly 


£  (  0  cV\ 

de  \tX’  0 ) 


from  (1.6),  so  (5.2)  becomes 


(5.1) 


(5-2) 


(5.3) 


i.c. 

— =  0  (^•^) 
The  equilibrium  condition  #  s=  0  becomes  IiHi  —  cX/tj  =  /aMa  ~  or  equivalently, 

U>1  as  U>3. 

Thus,  the  equilibria  are  given  by 

(i)  Ct>i  s  (<;2  0  or 

(ii)  wi  B  wj  ^  0,  X'  s  0 

Let  us,  for  simplicity,  choose  our  reference  configuration  so  that  dj2  ^21  are  parallel.  Then 


X'(f)  B  di2*d2isinf 


so  the  equilibria  in  case  (ii)  occur  when 
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(ii)'  either  d^a  =  0  or  dai  “  0, 
ot  9  —  0  or  IT. 

The  case  9  =  n  corresponds  to  the  case  of  folded  bodies,  while  9  =  0  corresponds  to 
extended  bodies. 

The  first  step  in  the  energy-Casimir  method  is  to  realize  the  equilibria  as  critical  points 
of  H  +  C,  where  H  is  given  by  (l.7b)  and  C  =  +  fti). 

One  calculates  from  (5.2)  and  (1.7)  that 


€yuiUJ2 

dE 


where 


(:)  =  = 


_1 //aMx 


The  first  variation  is 


<I(H  +  C)  =  +  (1^  +  *■)  4  (1^  +  ♦■)  (5.6) 

from  which  it  s  clear  that  critical  points  of  E  C  correspond  to  equilibria  of  (4.5)  with 

■  -(S.  -  -(S). 

where  the  subscript  ‘e’  means  evaluation  at  the  equilibrium.  As  in  other  examples  (the  rigid 
body  and  heavy  top  in  Holm  et  al.  [1085]),  ♦”(#*,)  is  arbitrary. 


The  matrix  of  the  second  variation  b 


<*(F  +  C) 


ft? 


§t9tn  (jia  _ 


I 


where 


e*H 

a*w 

5**Ta*ia 


a*w 

«Mr 


/2  - 
-£A  /x 


-  .Au,.) 

S^Jff  €»',  ,  ,  , 

3i^-  -^(-‘»"j  +  ^i“i) 


And 


a^H 


— eA 


,dn  dE' 

dm  dm. 


X.  „  d’jj 

=  -<Aa,.a„-.A^.,  -aa.,gjg- 

At  equilibrium,  A  =  ididj  (+  if  =  0,  —  if  =  ir)  bo 

1 


J-»  = 
d^E 


h  Ttdxdj 

(/|/j-€»d|d|)[q:£dicf3  /i 


and 


_ _  ^  _  d^E 

dSdm  SBm 


-Q^  =  -cA"u;’  =  ±cdidzu>l 


where  u,  =  wi  =  uj  ^  0  t.i  equilibrium.  Thus  (5.8)  becomes  becomes 

■±«Adi<i2<*'?  0 


£\E-i-C)  = 


a-+*»[; ;]] 


(5.9) 


This  matrix  is  clearly  positive  definite  if  di  0,  dj  0  if  9  =  0  (+  sign)  and  ^"(mc)  >  0 
and  is  indefinite  for  any  choice  of  4"(mc)  if  #  »  v. 


Another  way  to  do  the  stability  analysis  is  to  use  the  reduced  Hamiltonian  on  T*S^  given 
by  equation  (4.4).  After  completing  squares,  E  will  have  the  form  of  kinetic  plus  potential 
energy  with  effective  potential  given  by 


V(#) 


2A 


+  /»  -  *«A) 


+ 


(/i-/»)V  ] 

4(/i  +  fj  +  2cA)J 


(5.10) 


Minima  of  V  are  then  the  stable  equilibria  while  maxima  are  unstable. 

For  three  or  more  bodies,  this  method  of  looking  for  minima  of  the  potential  will  not  work 
because  the  symplectlc  structures  on  the  symplectic  leaves  will  have  magnetic  terms. 

Theorem  2  :  The  dynamics  of  the  2  body  problem  is  completely  integrahle  and  contains  one 
stable  relative  equilibrium  solution  ($  =  0  ~  the  stretched  out  ease)  and  one  unstable  relative 
equilibrium  solution  fff  =  ir  -  the  folded  over  ease).  The  dynamics  contain  a  homoclinic  orbit, 
as  in  Figure  1. 


|6  Multibody  Problems 

We  have  proved  that  the  Hamiltonian  formulation  of  the  previous  sectipns  extends 
in  a  natural  way  to  systems  of  N  planar  rigid  bodies  connected  to  form  a  tree  structure 
(Figure  5).  Since  the  general  statement  of  this  result  requires  significant  additional 
notation  and  the  explicit  introduction  of  the  notion  of  nested  bodies,  we  limit  ourselves 
to  the  special  case  of  a  chain  of  N  bodies  (Figure  6). 


Theorem  S  ;  The  total  kinetic  energy  (  Hamiltonian  )  for  cn  open  chain  of  N 
planar  rigid  bodies  connected  together  by  hinge  joints  takes  the  form, 

H  =  (6.1) 

where  p  =  {pi,P2,  -  •  •  <<  the  momentum  vector  and  3  is  the  corresponding 

NxN  pseudo-inertia  matrix  which  is  a  function  of  the  set  of  relative  (  or  joint  ) 
angles  between  adjacent  bodies.  The  reduced  dynamics  takes  the  form  : 


Ml 

M: 


SB 

S*i,t 

SB  _  SB 
Stt,»  St2,l 


M. 


SB  _  SB 


SB 


i  —  ag  _  SB 
where,  is  the  joint  angle  between  body  i  +  1  and  body  i. 


for  i=J,..,N‘l 


The  associated  Poisson  structure  is  given  by, 

If  [K.  JL)  ^ 

<-1  dg^+i,.  \dPi  3m.+i/ 


(6.2) 


(6.3) 


This  is  proved  in  a  way  similar  to  the  two  body  case  (see  Sreenath,  Krishnaprasad, 
Marsden  [1087]). 

The  structure  of  equilibria  and  the  associated  stability  analysis  become  quite  com¬ 
plex  and  interesting  as  the  number  of  interconnected  bodies  increases.  A  mixture  of 
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Figure  7 


topological  and  geometric  methods  may  be  necessary  to  extract  useful  information  on 
the  phase  portraits. 

In  the  remainder  of  this  section,  we  illustrate  some  of  the  complexities  of  multibody 
problems  by  giving  an  analysis  of  the  equilibria  and  stability  for  a  system  of  three  planar 
rigid  bodies  connected  by  binge  joints  (  see  Figure  7). 


§6.1  Three-Body  Problem 


The  Hamiltonian  of  the  planar  three-body  problem  is  given  by  equation  (6.1)  with 
the  momentum  vector  fi  and  the  coefficient  of  inertia  matrix  J  being  defined  as  below: 


H  =  (m.MJ.Ms)’’ 


?i  Aia(^j.i)  Asi(92,i  -f  ®s,j) 

Aij(^j.i)  It  Aasl^s,*) 

+  ^3.j)  Aaj(^s,j)  7s 


where,  7*8  and  A*8  are  defined  later.  #3,1  and  ^3,3  are  the  relative  angles  between  body 
2  and  body  1,  and,  body  S  and  body  2,  respectively. 


The  dynamics  of  a  three-body  system  of  planar,  rigid  bodies  in  the  Hamiltonian 
setting  is  given  by:  ^ 


Ml 

= 

BB 

0*2.1 

Ms 

= 

BB 

0*2.1 

+  i 

Ms 

= 

BB 

Btta 

^3.1 

BB  _ 

Bit) 

BB 

Bitx 

^3,3 

BB  _ 

Bit) 

BB 

Bit) 

BB 


(6.5) 


Remark  :  The  sum  (  Mi  +  momentum  variables  is  a  constant. 

Me 

Remark  :  The  coefficient  of  inertias,  and  A,y ’s  are  given  by: 

li  = 

7a  = 

7,  = 

^33(^3.:)  = 

^3i(^a.i  +  ^5.j)  = 
ti:  = 

where  mj’s  and  /,’s  are  the  mass  and  inertia  respectively  of  body  i,  and  d,y’6  are  deSned 
as  in  Figure  7. 

§6.2  Three-Body  Problem  :  Equilibria 

Refer  to  Figure  7.  Let  the  centers  of  mass  of  the  bodies  Oi.Os  and  Oi  respectively, 
also  be  the  origins  of  the  local  frames  of  references.  The  O33  be  the  joint  between  body 
1  and  body  2,  and,  Oaa  be  the  joint  between  body  2  and  body  3  respectively.  The  local 
coordinate  system  for  body  1  is  chosen  such  that  the  x-axis  is  parallel  to  the  line  joining 
Oi  and  Oil .  Similarly  the  coordinate  systems  for  body  2  and  body  8  are  chosen  to  be 


(A  +  (<w  +  «si)  <  dij,di3  >], 

[h  +  <  d3i,d3i  >  +  cas  <  das, das  > 

+  €si  <  (das  —  d3i),(das  —  dai)  >] 

[/s  +  (cjs  +  *s0  <  dsa»dsa  >]» 
{fiaA(_d„.d„)(^a.i)  +  f3iA(d„-dj,,d„)(^2.i). 
(c33A(-d„.da,)(^S.3)  +  esiA(_d„,da,-da,)(^S.3)), 
«SlA(_d„.d„)(^3,I  +  ^s.l). 

— .Dd  I.y  =  1.2,3, 

iTij  +  fUj  + 

x.ycos(o)  -h  (xxy)  Bm(a) 


parallel  to  the  line  joining  O2  and  Ojj,  and,  the  line  joining  Os  and  Ojs  respectively. 
DeHne  the  vectors  di2,d2i,d2s,d3],  in  their  respective  local  coordinate  systems  to  be 
;  di2  =  (ci,0],  d2i  =  (— fti.O],  d2s  =  ds2  =  |— di,0]. 

The  equilibria  for  the  three-body  system  can  be  found  by  setting  the  dynamical 
equations  in  (6.5)  to  be  zero.  This  results  in  the  following  equations  : 


Bn  _  BH  _ 

Bt».i  ~  Bt,,»  -  ^ 


Sz.i  =  a;2  —  Wi  =  0 
=  ws  —  W2  =  0 

From  the  above  equations  it  can  be  seen  that 


=  u>2  =  Ws  =  Wo(constant.) 

The  system  angular  momentum  ft,,  and  the  Hamiltonian  H  are  given  by 

=  Wo  [53  Ti  +  2(Ai2(52,i)  +  5^28  +  Xsi(»,  ,)(^2.1  +  ^S.l)) 

Lisl 


^  =  rwj  53/i  +  2(Aj2(^J,i)  +  5[23(^3.j)  +  A8i(^2,1  +  ^3,s)) 

^  Li=l 


(6.10) 


It  is  a  consequence  of  Theorem  3  and  (6.6)  that, 


2  5^2,; 


^ /v-l  T-1  \ 


1  dJ 


2(mi  +  m3  +  ms) 


[Ai  sin  (02.1  +  ^J,z)  +  Bi  sin(^2,i)  +  01  cos(^2.i)] 


or,  for  the  ooa-degeoerate  case  (a;#  0), 

Ai  tin(934  +  •s^)  +  Bi  1^(92, 1)  +  Oi  cos (02, 1)  *  0 
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(6.11) 


'  a  **4 


where, 


=  mim^eidi 

=  (mj(hj  +  e,)  +  mjhi]  miCi 


=  mimsCiej 


(6.12) 

(6.13) 

(6.14) 


Similarly,  for  we  get. 


5^3,2  2  (mi  +  mj  +  mj) 


[i4isin(52,i  +  tfs,j)  +  B2sin(tf3,2)  +  Cscos(^3,j)]  =  0 


where. 


Bj  =  (’’'1(61  +  Cl)  +  mjCijmstii 

Cj  =  (mi  +  m2)msrfiea 


(6.15) 


(6.16) 

(6.17) 


We  aasemble  the  final  equilibrium  equations  from  equations  (6.11)  and  (6.15)  : 

Ai  sin(53,i  +  ®5,j)  +  Bi  sin(93,i)  +  Ci  cos(fi2,i)  =  0 
j4i  Bin(^2.i  +  ^3.2)  +  Bi  sm(^3,2)  +  C2  005(^3,2)  =  0 


§6.3  Thrcc-Body  System  :  Special  Kinematic  Case 

We  consider  here  a  case  of  the  three- body  system  with  a  special  kinematic  structure 
where  the  centers  of  mass  of  the  bodies  are  aligned  with  the  joints  in  a  straight  line 
when  the  bodies  are  in  a  stretched  out  position.  In  this  case  we  shall  prove  that  (6.18) 
4  or  6  solutions.  For  this  situation  e  =  (cijej)**  =  {ei,0]^,  and  so  from  (6.14)  and 
(6.17) 


Cj  =  0 


Cl  =  Cj  =  0 


Thus  (6.18)  reduces  to 


Ai  sin(^2.i  +  ^s.j)  + 

Ai  sin(fl2,i  +  +  B2  sin(®s,j) 


(6.19) 

(6.20) 


Ai  =  cidimims, 

Bi  -  ((6i  +  ei)ms  +  himalcimi, 
B2  =  [(61  +  ei)mi  +  Cima]  dims. 


Subtracting  (6.19)  and  (6.20)  we  get 


where 


sm(^s.j)  =  Ksm(ff2,i), 


Bi 

““ft- 


Expanding  (6.19)  and  substituting  (6.24),  we  get 


i4i  sin(tfa.i)  lcos(tfs.j)  +  'ccos(tf2.i)  +  t)  = 


where 


(6.21) 

(6.22) 

(6.23) 

(6.24) 

(6.25) 

(6.26) 

(6.27) 


Consequently  from  (6.24)  and  (6.26)  we  have 


•in(tfa.i)  ®  0  and  ■in(ds,i)  =  0, 


(6.28) 


Figure  8  :  Fundamental  equilibria 


or, 

(6.29) 

cos{ff8.a)  +  »tcos(tfj,i)  +  T  =  0.  (6.30) 

It  is  obvious  from  coosidering  (6.28)  that  the  following  four  roots  of  the  pair 

could  be  readily  identified  : 

{0,0} 

.  (6.31) 

{ir.O) 

We  label  these  equilibria  as  the  fundamentat  equilibria.  A  stick  figure  representation 
(Figure  8)  helps  in  bringing  out  the  symmetrical  way  in  which  these  equilibria  occur. 

The  remaining  equilibria  for  this  system  are  computed  as  the  solutions  to  (6.29) 
and  (6.30).  Since  the  equilibrium  equations  are  nonlinear  and  parameter  dependent, 
one  needs  to  excercise  care  while  solving  them.  The  parameter  dependence  of  the 


equilibrium  solutions  can  be  summarized  by  two  sets  of  constraints  -  paramcter^sign 
and  parameter-value  constraints  respectively.  It  was  found  that  two  extro  equilibria 
(other  than  the  fundamental  equilibria)  can  exist  at  a  time,  subject  to  the  existence 
of  suitable  values  of  k  and  r  satisfying  these  constraints.  The  maximum  number  of 
equilibria  for  a  general  three-body  system  (special  kinematic  case)  is  thus,  6.  For  some 
values  of  k  and  r  not  satisfying  these  constraints  and  for  the  cases  with  k  and/or  r 
being  zero  these  extra  equilibria  merge  with  the  fundamental  equilibria  to  give  a  total 
of  4  equilibria. 

6.3.1  Parameter-Sign  Constraints 

This  constraint  set  restricts  the  existence  of  values  of  P&ir  depending  on 

the  signs  of  k  and  r. 

Using  (S.27)  in  (6.19)  we  get 

8in(tfj.i  +  ^s.j)  *  -r8in(^j,i).  (6.32) 

Taking  into  account  the  signs  of  k  and  r,  from  (6.29)  and  (6.32)  we  get  Figure  9,  which 
illustrates  the  feasible  regions  of  the  solution  pair  ^s,z)  to  form  the  parameter-sign 
constraints. 


6.3.2  Parameter-Value  Constraints 


The  existence  of  solutions  of  (6.29)  and  (6.30)  is  also  dependent  on  the  actual  values 
of  K  and  r  (which  are  constants  for  a  given  three-body  system).  The  parameter-value 
dependence  of  the  solutions  can  be  formulated  by  squaring  and  adding  (6.29)  and 


(6.30),  and  simplifying  to  get 


cos(*„)  =  -■  (6,J3) 

€Ot(*,j)  »  - — ^ - ,  (6.S<) 

•O. 

(6.35) 

(6.36) 


J 


K>0  t  >0 


k:<0  t  <0 


□  K>0  X <0 


K  <0  t  >0 


Figure  9  :  Parameter- Sign  Constraints 


These  equations  could  be  represented  in  the  form  of  a  graph  as  in  Figure  10.  The  graph 
has  been  drawn  for  k!  >  0  and  r'  >  0,  where 


k!  =  |k|, 
r*  =  |r|. 


6.3.3  Local  Frames  of  Reference 

It  is  necessary  to  choose  a  local  frame  of  reference  for  each  of  the  bodies  in  order  to 
parameterize  the  system  and  study  the  system  equilibria.  Refer  to  Figure  11.  Proper 
choice  of  the  local  frames  of  reference  for  bodies  1  and  3  results  in  the  vectors  c  ~ 
d  =  (d|,0]^,  where  both  Ci  and  di  are  positive.  In  general,  the  local  frame 
of  reference  of  body  2  could  be  chosen  In  such  a  way  that  c  —  (cii  —  [ei,0]^,  where 
Cl  is  positive.  Note  that  if  b  «  the  kinematic  parameter  bi  could  be  either 

negative  or  positive.  The  two  cases  of  the  signs  of  hi  represent  whether  the  center  of 
mass  of  body  2  is  inside  the  line  segment  joining  the  hinges  On  and  0%%  or  outside  it.  If 
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negative 


Figure  11  :  Reference  configuration 


aoy  of  the  kinematic  parameters  Cj  or  di  is  equal  to  sero  then  the  three  body  problem 
decomposes  into  a  two-body  problem  and  a  one-body  problem.  It  is  also  important  to 
observe  that  with  this  choice  of  local  frames  of  reference,  A\  is  positive  (see  (6.21)). 

6.3.4  Parameter-Dependent  Equilibria 

We  DOW  delve  into  particular  cases  of  the  signs  of  parameters  k  and  r  and  establish 
the  solutions  to  the  equilibrium  equations.  We  constantly  refer  to  (6.21)-(6.27)  while 
formulating  the  necessary  conditions. 

In  all  the  cases  we  consider,  we  first  ascertain  that  there  exist  physically  realizable 
values  of  the  kinematic  parameters  -  Ci,  6i,  ti  and  dx,  before  finding  the  actual 
solutions.  The  equilibria  are  evaluated  based  on  the  signs  of  cos(02,i)i  rind,  cos(0a,]) 
(see  (6.33)  and  (6.34)  ),  and  according  to  the  parameter-sign  and  parameter-value 
constraints.  The  results  are  presented  in  the  form  of  a  table  for  each  case.  The  graphs 
under  the  column  parameter-sign  constraints  have  to  be  read  with  ^3,1  as  the  X-axis  and 
4s, 3  being  the  Y-axis  (see  Figure  5.6  for  more  details).  The  shaded  regions  represent  the 
valid  regions  of  existence  of  the  {43,i,4s,3}  pair.  In  the  column  of  the  parameter-value 
constraints,  the  regions  referred  to  are  the  regions  of  Figure  5.7. 

Given  values  of  k  and  r,  one  can  identify  the  corresponding  table  depending  on  the 
signs  of  these  parameters,  and  determine  which  region  they  belong  to  with  regard  to 
Figure  5.7.  The  two  eztra  equiltbria,  if  any,  could  then  be  read  off  from  the  table. 

CASE  1 :  K>0r  >0 

Tot  k  and  r  to  be  greater  than  sero,  Ai,Bi  and  B3  should  be  greater  than  sero.By 
choice  of  the  local  frames  of  reference  we  have  from  (6.22)  and  (6.23)  : 

(hi  +  Ci)ms  ■4-  hxni|  >0  so  Ci  >  — 

(6i  +  ei)mi  +  cxms  >0  so  ci  >  -  ( 

\nii  +  ms/ 

i.c.. 


(6.37) 


Table  1  :  k  >  0  t  >  0 

This  is  automatically  satisfied  if  bi  >  0. 

The  equilibrium  solutions  are  given  in  a  compact  form  in  Table  1. 

CASE  2:  K<0.T  <0 

The  case  k  <  0  and  r  <  0  can  be  realized  if  and  only  if  Bj  <  0,  and  B]  >  0  (since 
>  0  always). 

Simplifying  so  from  (6.22)  and  (6.23)  we  have 

Naturally,  the  above  equation  indicates  that  this  case  is  possible  only  if  bi  is  negative 
(since  ei  >  0). 

Table  2  gives  the  equilibria  associated  with  this  case  if  (6.38)  is  satisfied. 


CASE  3:  K>0.r<0 


Cas0  ICospjj 


2.1  >  0 


2  2 


<0 


)  Cos{0j^ 


Parameter-angle 

oonstrainit 


Parameter  -value 
conttrarnts 


region  3 


region  1 


Equilibria 


2.3 


<0 


>0 


2  4 


>0 


<  0 


region  2 


not  satisfied 


Table  2  :  «<  0  t  <0 


For  this  case  since  >  0  we  have  to  have  Bi,  Bj  <  0  i.e., 


<  -(■«=)*■• 


With  the  choice  of  local  fr2Lnies  of  reference,  ej  >  0  so  this  case  is  possible  only  if  6^  is 
negative  and 


«i  < 


(6.39) 


The  equilibria  are  as  given  in  Table  3. 


CASE4 :  k<0.t>Q 

The  necessary  condition  for  this  case  is 


(e-ro) 


Case 


CostBjj)  I  Cospj^ 


Parameter -sign 
constraints 


Parameter  -value 
constraints 


Table  3  :  k  >0  t  <  0 

But  ej  >  0,  80  6i  has  to  be  negative.  Then  from  (6.40)  is  greater  than  1  but  less 
than  a  fraction  -  which  is  impossible. 

So  kinematic  parameters  satisfying  ic  <  0  and  r  >  0  can  never  exist. 
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